UNCLASSIFIED 
AD  NUMBER 


AD017205 

CLASSIFICATION  CHANGES 

TO: 

unclassified 

FROM: 

confidential 

LIMITATION  CHANGES 

TO: 

Approved  for  public  release,  distribution 
unlimited 


FROM: 

Distribution  authorized  to  DoD  only; 
Administrative/Operational  Use;  FEB  1953. 
Other  requests  shall  be  referred  to  Office 
of  Naval  Research,  One  Liberty  Center, 
Suite  1425,  875  North  Randolph  Street, 
Arlington,  VA  22203-1995.  Pre-dates  formal 
DoD  distribution  statements.  Treat  as  DoD 

AUTHORITY 

31  Dec  1965,  DoDD  5200.10;  ONR  ltr  dtd  26 
Oct  1972 


THIS  PAGE  IS  UNCLASSIFIED 


UNCLASSIFIED 


AD  NUMBER 

AD017205 

CLASSIFICATION  CHANGES 

TO 

confidential 

FROM 

restricted 

AUTHORITY 

ONR  ltr 

dtd  11  Dec  1953 

THIS  PAGE  IS  UNCLASSIFIED 


UNCLASSIFIED 


RepAaduced 
luf  the- 


ARMED  SERVICES  TECHNICAL  INFORMATION  AGENCY 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 


DOWNGRADED  AT  3 YEAR  INTERVALS: 
DECLASSIFIED  AFTER  12  YEARS 
DOD  DIR  5200.10 


UNCLASSIFIED 


2rr  _ 


CLASSIFICATION  CHANGED 


FROM  JRESTRICTED 


1 


7 


Insert  Class 


jCONFIDENTIALI 

Insert  Class 


1712 


°8|  . 


ScrrCHKR  1955 


By  authority  of 


-ainlMin  accordance  wtO* 

la  MM  tort  S Howmbw[ 

MUMbOptN’  . ; - 


This  action  was  rendered  by 


Name  in  full 


Date 


Document  Service  Center,  A S T I A 


RESTRICTED 

OFFICE  OF  NAVAL  RESEARCH 
CONTRACT  N6onr-27015 

ii 

THE  EFFECTS  OF  ROTOR  BLADE  FLEXIBILITY 

AND  UNBALANCE 
on 

HELICOPTER  HOVERING  STABILITY  & CONTROL 


By 

LEONARD  GOLAND 


®ulS  dot^ent  has  been  reviewed  in  accordance  with 
©FJfAVlNST  5510  ^reph.  5.  Kie  security 

/ciaafsif  ication  assiu.ru.  4.  is  corrwjt. 


£..t«  £ 


By  difefetion  of 
Chief  of  Naval  Resea 


PRINCETON  UNIVERSITY 

AERONAUTICAL  ENGINEERING  LABORATORY 

Report  No.  218 


RESTRICTED 


RESTRICTED 


Department  of  The  Navj 


Office  of  Naval  Research 


T£E  EFFECTS  OF  ROTOR  BLADE  FLEXIBILITY  AND  UNBALANCE 
ON  HELICOPTER  HOVERING  STABILITY  AND  CONTROL 


February  - 1953 


Prepared  by: 


LEONARD  GOLAND 


Approved  by: 


/Vo ! 1 

A. A.  NIKOLSKY 


RESTRICTED 


RESTRICTED 


TABLE  OF  CONTENTS 


1. 

Summary 

1 

2. 

Introduction 

2 

A 

3* 

Results  of  Analysis 

6 

V \ 

4. 

List  of  Symbols 

10 

5- 

Theoretical  Analysis 

13 

6. 

Procedure  and  Application  of  Method 

4i 

7- 

Effect  of  Blade  Flexibility  on  Automatic  Control  System 
Requirements . 

46 

0 

8. 

Figures 

49 

9- 

References 

57 

RESTRICTED 


1 . SUMMARY 

A method  is  developed  to  determine  the  control  response  of  a single- 
rotor helicopter  equipped  with  unbalanced,  flexible  blades.  The  detailed 
analysis  is  simplified,  vithout  great  loss  in  accuracy,  in  order  to  present 
a practical  method  for  use  by  the  helicopter  engineer.  The  detailed  analysis 
Is  included  for  those  interested. 

By  use  of  this  method,  the  influence  of  rotor  blade  flexibility  and 
unbalance  on  the  stability  and  control  of  a conventional  aircraft  is  investi- 
gated. Typical  values  are  assumed  and  their  effects  evaluated.  It  is  shown 
that  blade  flexibility  and  imbalance  can  have  a noticeable  influence  on  the 
damping  in  pitch  of  a helicopter.  Consequently,  the  period  and  divergence 
of  the  unstable  fuselage  oscillation,  as  well  as  the  automatic  control 
system  requirements  for  stability,  sure  affected.  Neglect  of  blade 
flexibility  and  unbalance  could  materially  reduce  the  effectiveness  of  an 
automatic  control  system  design. 
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2.  INTRODUCTION 


To  date , most  treatises  on  helicopter  stability  and  control  have 
dealt  with  the  problem  under  an  assumption  that  the  rotor  blades  are 
structurally  rigid.  This  was  necessitated,  principally,  for  the  sake 
of  simplicity  in  the  formulation  of  a basic  theory..  For  contrary  to  this 
assumption,  it  is  well  known  that  in  practice,  helicopter  blades  have  a 
high  degree  of  flexibility  and  are  subjected  to  deflections  by  the  actions 
of  various  inertia  and  aerodynamic  forces  encountered  in  flight.  Basic 
theories  have  now  been  established  and  refined  (i.e.  Ref.  1,  2,  and  3)5 
and  are  commonly  accepted  as  capable  of  predicting  the  nature  of  helicopter 
stability  and  control.  Therefore,  it  appears  appropriate  at  this  time, 
to  extend  the  basic  theory  in  order  to  study  the  effects  of  various 
existing  conditions  such  as  rotor  blade  flexibility  and  unbalance. 

The  motion  of  a conventional  articulated  or  "see -saw"  rotor  blade, 
with  reference  to  the  rotor  shaft,  may  be  described  as: 

(1)  Rotation  about  the  rotor  axis 

(2)  Oscillation  about  a flapping  hinge 

(3)  Oscillation  about  a blade  root  feathering  axis 

(k)  Structural  twist  about  a blade  span  axis  (elastic  axis) . 

This  axis  may  coincide  with  the  feathering  axis,  depending 
upon  the  structured  characteristics  of  the  rotor  system. 

Additional  components  of  the  motion  are,  oscillation  as  a pendulum  in  the  plane 

of  the  rotor  disk  about  a drag  hinge  (lagging),  and  blade  structural 

bending  in  the  flapping  plane.  These  additional  components  can  materially 
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effect  the  twisting  tendency  (unbalance)  of  the  blades.  However,  as  was 
pointed  out  in  Ref.  4,  experimental  evidence  has  shown  that  for  articulated 
blades,  these  components  directly  have  only  a second-order  influence  on  the 
air  forces  on  the  blade.  This  may  not  be  the  case  for  rotors  in  which  the 
blades  are  rigidly  connected  to  the  shaft,  or  have  blades  of  extremely  low 
bending  stiffness.  In  these  instances,  structural  deformations  (bending)  in 
the  flapping  plane  may  be  of  prime  importance  (Ref.  6).  The  analysis 
presented  herein  is  principally  concerned  with  helicopters  having  an 
articulated  or  "see-saw"  type  rotor. 

On  the  other  hand,  it  has  been  shown  by  Wheatley  (Ref.  4 and  5),  that 
the  torsional  flexibility  and  unbalance  of  a rotor  blade  has  a pronounced 
influence  on  the  flapping  motion  of  a blade.  Except  in  the  case  of  balanced 
blades  (i.e.  blades  having  symmetric  airfoils  whose  elastic,  aerodynamic, 
and  mass  gravity  axes  coincide),  the  resultant  of  the  air  forces  and  mass 
reactions  produce  a couple  which  tends  to  twist  the  blade.  This  resulting 
twist  has  a constant  magnitude  for  pure  hovering  flight,  while  in  disturbed 
or  forward  flight  it  has  periodic  values  as  well.  This  twist,  naturally, 
affects  the  flapping  motion  of  the  blade,  and  consequently  the  characteristics 
of  the  helicopter's  motion. 

The  existence  of  the  above  phenomenon  has  long  been  appreciated  by 
blade  designers,  due  to  its  associated  disturbing  effect  on  stick  forces. 

To  eliminate  blade  twist  and  the  disturbing  stick  forces  they  have  striven 
to  design  balanced  blades.  With  this  condition  satisfied  the  disturbing 
stick  forces  due  to  the  blades,  as  well  as  blade  twist  would  be  eliminated. 
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For  it  is  well  known  that  in  the  conventional  rotor-control  system,  oscillating 
stick  forces  are  due  entirely  to  blade  unbalance,  having  both  mass  and  aero- 
dynamic origins.  Unfortunately,  however,  the  attainment  of  balanced  blades  | 

in  practice  is  extremely  difficult.  This  fact  is  readily  apparent  when  one 
considers,  as  previously  stated,  blade  distortions  in  flight,  in  the  flapping  , ' 

and  lag  planes.  In  many  cases  a balanced  blade  design  proves  impractical. 

In  addition,  constructional  inrperfections  and  use  often  cause  unbalance. 

Consequently,  blade  twist  and  oscillating  stick  forces  axe  present  in 
varying  degrees  in  helicopters. 

Due  to  the  disturbing  nature  of  these  oscillating  stick  forces  recent 
helicopter  designs  resort  to  the  use  of  irreversible  control  systems.  This 
approach  tends  to  eliminate  the  disturbing  stick  forces.  However,  since 
the  cause,  blade  unbalance,  has  not  been  eliminated,  blade  twist,  and  its 
associated  effects  are  still  present  in  these  designs. 

i 

Contrary  to  the  above  tendency  to  eliminate  automatic  blade  twist, 
there  is  a feeling  among  many  that  effective  stability  and  control  may  be 
obtained  by  proper  design  of  the  rotor  system  utilizing  unbalanced  blades. 

Naturally,  such  a system  would  entail  a minimum  weight  penalty  on  the 
already  relatively  poor  carrying  capacity  of  helicopters.  For  unlike  the 
fixed-wing  aircraft  where  automatic  stability  and  control  requires  the 
installation  of  special  gyroscopes,  the  helicopter  can  make  use  of  its 

I 

i 

rotor  system  which  itself  is  a gyroscope.  \ 

1 

i 

j 

The  results  which  follow  attempt  to  give  some  indication  as  to  the 

o 

effects  of  blade  flexibility  and  unbalance  on  the  stability  and  control  i 

t 

! 
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a typical  helicopter.  The  effects  on  the  automatic  control  system  re- 
quirements for  stability  are  discussed  in  Section  7. 


RESTRICTED 


5 


RESTRICTED 


3.  Results  of  Analysis 

Similar  to  fixed-wing  aeroelastic  phenomena,  the  effect  of  blade 
flexibility  depends  primarily  upon  the  relationship  in  flight,  between  the 
blade's  elastic,  aerodynamic  and  mass  gravity  axes.  Using  the  analysis 
presented  in  Section  5,  both  mass  and  aerodynamically  unbalanced  blades 
have  been  investigated.  To  demonstrate  the  procedure  and  application  of 
the  analysis,  a typical  calculation  is  presented  in  Section  6. 

For  the  quantitative  investigation,  the  parameters  of  a typical 
5000  pound  helicopter  were  assumed  (section  6)  and  the  results  compared 
to  that  obtained  under  the  assumption  of  structurally  rigid  blades.  The 
blades  were  considered  fixed  in  torsion  at  their  root  end,  and  to  be  of 
uniform  spanvlse  construction,  thus  having  principal  torsion  modes  of  the 
typo  S (r)  ■ sin  #7*.  It  is  worthy  to  mention  however,  that  in  the  analysis 

W 

(Section  5)  >>1000  twist  is  introduced  through  a general  torsion  mode,  # (r), 
and  consequently  the  analysis  may  be  applied  to  helicopters  with  blades 
having  any  arbitrary  torsional  stiffness  or  mass  distribution. 

a)  Mass  Unbalanced  Blades 

For  overbalanced  (center  of  gravity  forward  of  the  elastic  axis)  as 
well  as  underbalanced  (c.g.  rear  of  elastic  axis)  blades  the  controlling 
dimensionless  flexibility  parameter  was  found  to  be  j ; where, 

12  is  the  blade  mass  product  of  inertia  about  the  flapping  hinge  and  blade 
elastic  axes  (positive  for  forward  c.g.)^  jfL  is  the  rotor  speed,  aadK  is 
the  blade  torsional  stiffness. 
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The  principal  effect  of  overbalanced  blades  is  to  increase  the  damping 


in  pitch  of  the  helicopter.  This  is  shown  in  Fig.  1,  in  which  the  Increase 

in  damping  in  pitch  of  the  helicopter  is  given  for  various  values  of 

relative  to  the  damping  for  = 0.  This  effect  can  be  explained 

as  follows.  With  conventional  rotor-cohtrol  systems,  the  rotor  tip-path  plane 
tends  to  pitch  with  the  helicopter  fuselage.  However,  in  so  doing,  coriolis 
forces  which  oppose  pitching  of  the  tip -path  plane  are  created  on  the  blades. 
These  forces  are  proportional  to  the  pitching  velocity  of  the  tip-path  plane, 
which  is  approximately  equal  to  the  fuselage  pitching  velocity.  Thus,  pitching 
of  the  tip-path  plane  lags  behind  the  fuselage  pitching  motion.  The  resulting 
inclination  of  the  tip -path  plane  relative  to  the  rotor  shaft  is  responsible 
for  the  pitch  damping  of  the  helicopter.  Moreover,  the  above-mentioned 
coriolis  forces  when  acting  on  unbalanced  blades  cause  the  blades  to  twist. 
Overbalanced  blades  twist  so  as  to  increase  the  inclination  of  the  tip -path 
plane  and  thereby  increase  the  pitch  damping.  This  effect  is  therefore  quite 
similar  to  that  obtained  by  increasing  the  "heaviness"  (I  ) of  the  blade. 

The  opposite  is  true  for  underbalanced  blades. 


The  response  of  the  helicopter  in  pitch  to  a control  step  input  is  shown 
in  Fig,  2 for  various  values  of  . The  response  in  the  first  second 

is  shown  more  accurately  in  Fig.  3.  It  is  noticed  that  as  the  blades  become 
more  flexible  or  more  overbalanced,  decreased  fuselage  oscillations  of  longer 


period  result.  This  comes  about  principally  through  the  increased  stability 
in  pitch  of  the  craft.  However,  as  pointed  out  in  Section  7,  without  control 
displacements  in  phase  with  fuselage  attitude,  the  helicopter  cannot  be 
stabilized  by  increased  pitch  damping.  In  the  limit,  only  a neutrally  stable 


craft  results.  The  damping  rates  and  frequencies  of  the  unstable  fuselage 


motions  of  Fig.  2 are  as  follows: 


Damp,  factor  (m^) 


Frequency 


Time  to  double  Period 
amplitude  (sec)  (sec) 


0 0.188 

.2  0.179 

1.4  0.070 

2.2  0.045 


0.441 

3-7 

14.2 

0.441 

3*9 

14.2 

0.378 

10.0 

16.6 

0.344 

15.6 

18.3 

In  the  case  of  underbalanced  blades,  increased  fuselage  oscillations  as  well 
as  decreased  periods  are  to  be  expected  due  to  the  decrease  in  pitch  damping. 
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By  the  addition  of  a proper  amount  of  damping  to  the  twisting  motions 
of  the  blades  it  is  possible  to  introduce  a time  lag,  and  the  blades  could 
effectively  have  a response  component  in  phase  with  attitude  as  well  as 
pitching  rate.  The  amount  of  torsional  (structural  and  aerodynamic) 
damping  introduced  by  most  blades  is,  however,  relatively  small. 

Blade  flexibility  and  unbalance  can  materially  affect  the  automatic 
control  system  requirements  for  stability.  In  fact,  consideration  of  blade 
flexibility  and  unbalance  could  mean  the  difference  between  satisfactory  or 
unsatisfactory  behavior  of  a automatic  control  system.  This  effect  is 
discussed  in  Section  7* 
by  Aerodynamlcally  Unbalanced  Blades 

In  the  case  of  aerodynamically  imbalanced  blades  the  controlling 

2 3 

dimensionless  parameter  was  found  to  be  Ch/K,  where  Ch  = 1 /2  (f  ac  12  R h^. 

"<f"  is  the  mass  density  of  air,  "a"  the  aerodynamic  lift  coefficient  of 
the  blade  and  "h^",  the  aerodynamic  chordwise  unbalance  in  feet  (positive 
for  forward  of  elastic  axis).  The  other  symbols  are  as  previously  defined. 

Similar  to  the  action  of  mass  overbalanced  blades,  aerodynamically 
underbalanced  blades  increase  the  damping  in  pitch  of  the  helicopter. 

Fig.  4 shows  the  relative  increase  in  pitch  damping  of  the  helicopter  for 
the  case  of  aerodynamically  underbalanced  blades.  This  increase  can  be 
attributed  to  the  aeroelastic  twisting  of  the  blades.  Similar  to  the  action 
of  mass -overbalanced  blades,  aerodynamically  underbalaneed  blades  twist  so 
as  to  increase  the  inclination  of  the  rotor  tip -path  plane  relative  to  the 
pitching  rotor  shaft..  Aerodynamically  overbalanced  blades,  on  the  other 

hand,  decrease  “the  inclination  and  consequently  the  pitch  damping  of  the 
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The  response  of  the  fuselage  in  pitch  to  a control  step  input  is  shown 
in  Fig.  5 for  various  values  • The  response  in  the  first  second  is 

shown  more  accurately  in  Fig.  6.  As  would  be  expected  with  increased  pitch 
damping,  the  amplification  of  the  fuselage  oscillation  is  decreased  and  the 
period  is  increased.  The  damping  factors  and  frequencies  of  the  unstable 
oscillations  of  Fig.  5 are, 


Ch 

K 

Damp,  factor 

(n»l) 

Frequency^) 

Time  to  double 
amplitude  (sec) 

Period 

(sec) 

0 

0.188 

o.44i 

3-7 

14.2 

-2.1 

O.I36 

0.422 

5.2 

14.9 

"3-5 

0.106 

O.398 

6.6 

15.8 

For  the  uniform  blade  considered,  the  amount  of  unbalance  represented  by 
the  above  values  of  Ch/K  are, 

i Ch/K  aerodynamic  unbalance  (#  chord) 

0 0 

-2.1  5-5* 

-3-5  9.1# 

As  in  the  case  of  mass  unbalance,  aerodynamic  unbalance  blades  do  not 

respond,  relative  to  the  rotor  shaft,  to  fuselage  attitude  but  to  the  rate 

of  pitch  of  the  machine.  Consequently,  the  effect  on  the  automatic  control 

system  requirements  for  stability  is  similar  to  that  of  mass  imbalanced 

1 

blades  (Section  7) • 
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LIST  OF  SYMBOLS 

l ' rn 

a = aerodynamic  lift  coefficient  of  the  "blade ^ 

2 

A » rotor  disk  area,  ft  . 

&2  m l/2pbc.TlR2 
b = number  of  blades 

B1  = control  angle  of  attack  input  to  blades  (sine  f component),  radians (Ea. 4) 
c = blade  chord,  feet 

Ch  = blade  aerodynamic  unbalance  parameter,  ft.  •//  — ~TL 

E,  L,  N,  etc  = blade  motion  influence  factor*  (see  Eq.  59) 
g ■ acceleration  due  to  gravity,  ft/sec^ 

h = distance  of  rotor  above  helicopter  center  of  gravity( including  blades), ft. 

h^  * aerodynamic  chordwise  unbalance,  positive  for  forward  a.c.,  ft. (Fig. 8) 

Iy  = moment  of  inertia  of  helicopter  about  a lateral  axis  located  at 
helicopter  c.g.  (including  blades),  the  mass  of  the  blades  being 
considered  concentrated  at  the  rotor  hub,  slug  - ft^ 

K * blade  torsional  stiffness  (ft. ///radian) . 

UT  = mass  of  helicopter,  including  blades,  slugs  *■*  0 

R = effective  aerodynamic  blade  radius,  ft. 

r * radial  distance  from  rotor  hub  to  a blade  element,  feet  (Fig.  7) 
t - time,  sec. 

T = rotor  thrust,  pounds  — • >/ 

W * gross  weight  of  helicopter,  pounds 
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fl  » rotor  angular  velocity,  radians  per  sec. 

= blade  azimuth  position,  radians  (Fig.  7)  < 

0 = flapping  of  blade  neutral  axis  from  horizontal  plane,  radians  (Fig.  7) 

* jS0  - a.,  cos  p - 6,  s/s?  j 

I 

o<,  = fuselage  pitch,  radians  (Fig.  7)  j 

JS  = chordwise  distance  of  blade  center  of  gravity  from  neutral  axis,  J 

positive  for  forward  c.g.,  ft.  (Fig.  8) 

^5(r)  « blade  torsional  mode  shape  (see  Eq.  4). 

0(r)  = geometric  angle  of  attack  of  a blade  airfoil  section  (Fig.  8)  ! 

| 

©Q  = cpntrol  collective  pitch  of  the  blades 

to  a blade  twist  at  tip,  collective  component,  radians  (Eq.  4) 
t,  a blade  twist  at  tip,  cosine  JP  component,  radians  (Eq.  4) 
jQ  a blade  twist  at  tip,  sine  <0  component,  radians  (Eq.  4) 

* mass  density  of  air  ( slugs /ft3) 

i 

Ac  a dimensionless  rotor  induced  dovnwash  j 

F > dimensionless  blade  parameter 

dimensionless  horizontal  linear  velocity  of  rotor  hub 
* 

a Xp 

SL  R 

A « complex  roots  of  characteristic  equation  for  fuselage  oscillation, 
l/sec. 
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Rotor  Blade  Inertia  Integrals : 
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5.  THEORETICAL  ANALYSIS  OF  HELICOPTER  HOVERING  STABILITY  AND  CONTROL  INCLUDING 
THE  EFFECT  OF  ROTOR  BLADE  FLEXIBILITY  AND  UNBALANCE 

Method  of  Attack. 

The  general  method  of  attack  is  quite  similar  to  that  presented  by  Nikolsky 
in  Ref.  1,  in  which  however,  the  blades  were  assumed  to  be  structurally  rigid. 

In  the  above  reference,  the  helicopter  system  was  considered  to  possess  three 
degrees  of  freedom,  namely,  fuselage  horizontal  translation  and  pitch,  and  blade 
flapping.  In  this  analysis  to  account  for  blade  twist,  an  additional  degree  of 
freedom  is  introduced.  Blade  twist  is  introduced  oy  assuming  a general  blade 
torsion  mode,  (f)  (r) . This  permits  the  results  of  the  analysis  to  be  applied 
to  helicopters  with  blades  having  any  arbitrary  torsional  stiffness  or  mass 
distribution.  Basically,  four  equations  of  motion  are  developed  to  describe  the 
resulting  behavior  of  the  helicopter.  To  preserve  the  continuity  of  this 
development  certain  parts  of  Ref.  1 are  repeated  herein. 


Although  in  the  following  investigations  the  motions  axe  assumed  to  take 
place  in  the  longitudinal  plane  of  symmetry  of  the  helicopter,  "the  analysis 
equally  holds  for  the  lateral  motions.  For  simplicity,  and  since  this  analysis 
is  principally  concerned  with  the  effect  of  blade  flexibility  and  imbalance, 
the  case  of  a helicopter  with  no  effective  flapping  hinge  offset  is  treated. 

However,  by  use  of  Ref.  1,  and  the  following  procedure,  the  effect  of  a 
. flapping  hinge  offset  may  be  included.  In  addition,  the  tip  loss  factor  is 
i accounted  for  herein,  by  using  an  effective  aerodynamic  rotor  radius  (R)  in 

all  terms  except  those  involving  inertia  effects,  such  as  the  blade  inertia 
integrals  (Lp  Iq,  etc.,  pg.12  ). 

1 
! 
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Coordinate  Axes 

The  coordinate  systems  are  shown  in  Fig.  7 and  8.  The  stationary  vertical 
axis  is  assumed  to  coincide  with  the  rotor  shaft  when  the  helicopter  is  in 
its  undisturbed  hovering  condition.  The  stationery  XQand  Y0sxes  are  perpen- 
dicular to  the  vertical  axis  to  form  a right-handed  set  of  axes. 

The  moving  blade  axes  (x  and  y)  are  taken  to  coincide,  respectively, 
with  the  blade's  neutral  and  flapping  axes. 

Allowing  Xjj,  y^,  and  Z ^ to  represent  the  instantaneous  coordinates  of 
a blade  element  with  reference  to  the  stationary  axes,  the  coordinates  and 
their  time  derivatives  are  as  follows, 

*/>  - COS  fa  y-  JF  S'//?  fa 

fa  m - cs/s?  fa  - jZ cos  fa 

z*  - fir  v- 

Xa  — Xo  + rJ~L  s/s?  fa  y-  jfSL  cosy; 

fa  - cos  fa  s/s7  fa 

• - Ba  * fir  y- 

~ jfo  r r-Ihf  cos-  fa  — g sr/s?  fa 

fa  ~ r -flfi  s/,^  fa  g cos  fa 

• ••  •!  _ •• 

Z6  - fi>/~  -/■  g & 
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where  the  dots  refer  to  the  time  derivative  of  the  function. 


The  equation  of  motion  of  a blade  about  its  flapping  hinge  is  obtained 
by  equating  the  total  moment  at  the  flapping  hinge  to  zero.  Considering 
moments  positive  that  tend  to  increase/^  , this  may  be  stated  as. 

My  = 0 = (My)m  + (My)a  (2) 


where 

(My)m  = the  moment  due  to  the  inertia  loads. 

(My)a  = the  moment  due  to  the  airloads. 

The  moment  due  to  the  inertia  loads  may  be  evaluated  from  the  expression, 

ft 


(My)u  = -J [escos'/Ss-  v-  ca?  ¥■ 

-O  ) (3) 

* y 's/^  <£’G>)ja''r7 

Following  the  usual  convention,  the  blade  flapping  angle  and  blade  angle  of 
attack  may  be  expressed  as, 

& -a,cas  p — <4  jsy/v ^ 

~<Cco3  ~cT2  4? 


where 


cQ  *■  &o  -*•  to 
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where  tft  (r)  is  the  principal  torsional  deflection  mode  for  the  blade, 
is  the  constant  value  twist  of  the  blade  tip,  while  )f,  and  are, 
respectively,  the  cosine  and  sine  components  of  the  blade  tip  twist.  From 
Eq.  (4),  the  following  flapping  and  blade  twisting  velocities  and  accelerations 
are  obtained: 


— - (a,  +£  </z  - fC  - ^ 

3 — — {cl,  Z CL  — a,  ^TL^Jcas 


~ ('^z  — ^ et./  -CL-  — -CLC J s-/s7 


- -(<C  -ClJccxs  <//  ~^4l  -J7 

— ~ ftC/  -z  L?  £CL  — cC -CL? J eras- fZ 

~(*£  - CcC-CL  - J2  -CL-ysys?  ft/ 


As  j9  is  relatively  small  for  conventional  aircraft,  cos  @ fit/  1,  and 
the  first  term  on  the  aright  hand  side  of  Eq.  (3)  may  be  expressed  as. 


' £Ejk  /'cZ/r?  **/  f0  /-*c/s?7  Z f J; eZ/T? 


(6) 


where 


J zc  4 -Cl  -cz,  jnfjc. os- (7) 

~-T,  ^4  - - 4 CL  £Js-//?  fZ 
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o 


~[^3  K A ShJ<,+£jr3_f2  <£  ^ 

-j rs  ji_*  rjcos  ^ -fe  4 

-£jr3s2.  r,  - — t3  _CL*<&J^s?fs 


The  second  terra  on  the  right  hand  side  of  Eq.  (3)  may  "be  expanded  to, 


. ~./r 

' X4C0S  y/S-pc/tr?  -/■  / COS  &c//r? 


(9) 


whose  values  are. 


f X*>  cos  </'/•/£  c//r?  «•  X,  co>5  ^^X  — JT£  jn?(&o-5~//7 ^ c.os^ 

° +X?  xi/60  ocsps  - %cX^  -Tl2 d.  / car  </>  -^'XJX7A/  s/s? ft* 

jm^tsc,  jr//? y^X^.  X2~  ^ c<osr^' 

jx2 cos^'Xf&o'/??  = J~#  xT* &0  cos  xJT3X2-£  tfocos*^ 

0 “-^P  <&  S-//7  y/cos  pX  7*  -Z>  j£o  <&>  cas^y-X^  J£cosf* 

~XZ,  X2?Q  <r#7f'cdsfS'  ~^-JZ~3  XX^tfcos  ^ 

-fr-T*  -r^Y^  x^js/s?  ^ - A?  -z;  ~n  *&  s//?fxL1) 

+ Jfr-  JTs  sys? px  Y"  j7?7Tp  XX* (7%  Cos ^ 


RESTRICTED 


17 


RESTRICTED 


Similarly,  the  final  term  on  the  right  hand  side  of  Eq,.  (3)  may  he  expressed  as, 


R, 


Yb  sin  ^ dm  + f y sin^/'^  dm 


(12) 


where, 
,/P 


jTlf- ceu.pt'sys?^  (13) 

2-  ^ -3^  ^ 

jr/s?  ** 


j&P&s?  (/Sc//??  ■»  _y2.^Sfe  7*  -2£ -JO.1* (£ ^ysr  *$£/ 

y-JT^  czasp^sy^p*  -+-Z&  -O-^  <K>  cc>x  (/'jr/s?  p// 

c&J  </*  — Zfi~ -Za  -+  ^Jsv*?  </> 

~dfT -Z^  -/^ <£_  jr/y*  pt  — z^r  JQ?  _S2.£  sys?  p& 
-^Jccsp  -zjrJ&-f2*  fecaspf 


1 


l 
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Thus,  the  moment  due  to  inertia  loads  becomes, 

bj  J*,) 

S2-  J ~jf  s*Z 

- yXL  fcjcc&fb  Y-Zy'^ -^ea.,^1.)  (15) 

j=Fp  ^ ^ -T2-  <F/J  J s/*?Y 

--f2.aY~  ~ =?  SI*  - —2.^  fo 

> -A-/  -A-/ 


The  moment  due  to  the  air  loads  is 


where, 

- U/=>  **-2-6  Co^/S  -A-(Oq  cojs  ^ ~ c\z.  - 

~-^o  X0Jcz>sfb/  *-  ('—b'-TZ-C*.,  7 

•/  ^«sc  -b~b/? 


C/r  = JTjt-S/S?  ys— y^C-CZ^r /; 

— ' Ac,  3-//?  </s  ■/■  /^.JZL 


and 
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Consequently, 


and 


£ 

Ur  = ft* 

@Sr  ~ f&o-*  <£T  S^Slf  (K  ccxs-ps 

(17) 

+fzA>f-n£Q>+  ^ y.  p 


Thus,  the  first  integral  of  Eq.  (16)  becomes, 

Jeuft-cfr  -+s%_n?£  -Sj-n?  r,c.osp 

+{dts*  u qs?*g±  +<?S£x»  Si  & 

- -/2*  djP*,  -S?S2f  tfejsy/? ^ 


which,  may  finally  be  written  as, 

-/p 

>» 


±£2£SS^_  _ Z££r 


■+(§-&?<&  +3: 

&)**»& 


>To  y* 

°°  ~J 


(18) 
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Evaluation  of  the  second  term  of  Eq.  (l 6)  leads  to, 

U/=>Ur  A*.  -/2.V/P 

yi‘  ^ ^ ^ yt-sAa-.  U2.  /F’U^y)  s/s? 


(-19) 


/r  . 

( Up  Ur  f'd'f  — A**-  d-~~  y. 


+ * 


(20) 


- dZSfjcca  p +f.y£?_s2_*-a,  £ 

>3  , 


"*  Aa-  — A J -sy*? 

&[*<.£- ^Up-  Ur-Z'OA  =T  /^  A.  V4  JIt  -Jvif^3=25 - ^ (21) 


Substituting  Eqs.  (15)  and  (l6)  into  Eq.  (2)  results  in  an  equation  whose  terms 
are  dependant  on  sin  Y } cos  or  are  independent  of  these  harmonic  functions . 
Therefore,  in  accordance  with  Fourier's  analysis  this  equation  gives  rise  to 
three  independent  equations,  namely,  constant  terms  (steady-state), 
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sine  terras, 


s' 


There  has,  consequently,  heen  obtained  three  equations  to  describe  the 
resulting  blade  flapping  motion. 


Blade  Torsion  Equation 

The  equation  of  motion  of  the  blade  about  its  neutral  axis  may  be 
developed  by  evaluating  the  moments  about  this  axis  due  to  inertia  loads, 
airloads  and  elastic  restraint.  Considering  moments  positive  that  tend  to 
decrease  the  angle  of  attack,  the  equation  may  be  stated  as, 

Mx  * 0 * (Mx)m  + (Mx)a  + (Mx)d  (2 

where, 

(Mx)^  = the  moment  due  to  the  inertia  loads 
(Mx)  = the  moment  due  to  the  airloads 

B. 

(Mx),  » the  moment  due  to  the  elasticity  of  the  blade, 
d 
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The  moment  due  to  the  inertia  loads  (Mx)m  may  he  evaluated  from  the 

expression: 

/r» 

S'/S?  sjsjf-  Xs  sv/7 

0 ..  . 7 (26) 

7 cos  ^Jf'&/c//s? 

Evaluating  the  various  terms  specified  above  there  is  obtained, 

J'x*,  Jf  e/?s?  XJf  <//S7  c//r) 

~ -G-  _OlJ cos  ('J,  JTL 

~ -/£?  JT2J? jcos  JZT-^c^y  s-Z^ 

-zj:a  jn-  -JT?-n?c*<,  -^z^.z^Js/s? 


4&  cos  ^jt^Z/??  = _Z"=,  -Til t&a  —_ZT7  -ZzfjZb  s/s? f/>  <sos 


jQa.,<?os  ^ -z 


—fZr  Z~7  <£2,  s'//?  ^ -/■  sf~-Z~-,  _f2  /£>/  craxs  ^ 


( c/+ s/s?  ^^/Sofr?  — -Z^.  Z2f/£o  s/s?  v SZZ,  Z~lf/^„  ~s/s7  cos  (// 

a,  cos p-j^zitz:4s/*f/  (29) 

~ '^T-X2.2'-Z7  d,  s/s?  ^ — ^-ZZ,  uT2-'z^},  CSS- 
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./f 


(30) 


pi  &<//*?  “ -5-^  ^ -s -//?/' -jr7  jzL£*=?a&x>y> 

~~-Z~^  -^J2  tf^o"  —^rZZZ^  _/2'2  ^ 5~^  ^ 

^ -S*-*  fa  cts*^  -s£-^ZL  uG^Z^ctsc  fas  -^~-4  ~GMfa  <r<ar fa 

~Gf  ZZZ?  _/2.‘2^  cress-  y-  fafa-Zfa  —fLZ^z?  sr/s?  fas 
ZZ'j,  J2?<=*,  s'//? -0~z  fa  s-/s?  fas 

locos' y/J&otv?  =r  -^ZZ£  far/s?  fas  cos  fa  (3D 

> 

“*^■-22  _nffa  <zas fa  ^_z^_/q 

-*-ZZ^  _J2?  facts***  fa  — facos fa  ~^?r-ZZ,  s//?fa 

—far-ZZg  ~/2~*  fa  jts/9  fa 

Combining  Eqs.  (27)  through  (31)  as  specified  in  Eq.  (26),  the  moment  due 
to  inertia  loads  becomes  equal  to, 

fokL  =zfa,-n?&o  +[-zfa(£^4_n.j 

-jr#  fa  -^?_ZV  X.^- _/2  ^ + jz.Te.a-  fa 

-/  JZTs.  fa  faejc^  fa  -'■Jr-2*  f4  — <£?&./  -G-J -_ZZ,  jfa  (32) 

-*i£,  --&£  <%*.  +zTefiXlJs>ny 
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The  moment  due  to  airloads  has  the_yalue, 

fab  =■  ^ £/rJ<S/~ 

O S 

This  may  he  evaluated  simply  hy  use  of  Eqs.  (17)  and  (19)  to  yield, 

■fad.  - j~ =aw^~&c,  •+ 


< 3*  j» 

-J&  + Au ?3:  /k2tt£.% 

/ 3* 

Au  JO.  /P 


(33) 


(34) 


The  moment  due  to  the  elasticity  of  the  hlade  may  he  written  as, 

faxje/  ~ Sr'f'/o  - ce>*  ^ ^ svs?  (35) 

where  K is  the  effective  torsional  stiffness  of  the  hlade  and  root  attach- 
ment combined. 


The  hlade  torsion  equation  may  now  he  determined  hy  substituting 
Eqs.  (32),  (33),  and  (35)  into  Eq.  (25).  Similar  to  the  case  of  the  flapping 
equation,  the  resulting  torsion  equation  contains  terms  which  are  dependant 
on  sin  p/ , cos  jp  , or  which  are  independent  of  these  harmonic  functions. 
Therefore,  this  equation  also  gives  rise  to  three  Independent  equations, 
namely,  constant  terms  (steady-state), 

/rvs  i-jr7^L2^  ra 


*JG>S 

&o 


(36) 


-/  ^ -J2-2  tz  y -r  O 
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and  since, 


/t/-  f>CLC^>-n^/r  ^ ~3r  A<^  v- 


Lx*  i 

~JEp3r  ooy 


yi;  -CD 


cosine  terms, 


“V^"  <^T  — JZIg  dj  —^JZZ, . 

y*  ^:-'lZ^  —fZL.  it-  JZ^£~^C>  S*o 


* 


-7' 

Z?, 


•4~.,^L*.  — < ***^0  ' j£o  7_ 

3>-n-  $5?'*  J 


where  Ch  = _ pac 
2 


JlW. 


(37) 


(38) 


(39) 


sine  terms, 

-A"  - J^(£,  - B<3.,-fl-J  -JT-,  c£y  -_Z^)  £ t LZJ^S2  <%  ^ 

~~^Z^r~&0/X0  <%X0  — _r^/r^ 

J'  o'  zf3*-?  ^ 5*  V -a./f’J 


Equation  of  Horizontal  Translation 

The  motion  of  the  helicopter  parallel  to  the  xQ  axis,  assuming  that®*; 
is  small,  may  he  described  by  the  equation, 

ir?  fx0  - /?c*J  -y-  M*a  = ° 

where, 

Hx  = the  con?) one nt  of  the  left  and  drag  forces,  acting  on  the 
o 

blades,  along  the  xQ  axis. 
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2.z>-  ^./p 


: ^ f3^  Or  cfs-//7  ^ - a.^  C4-Z+  LfaOrJ 

° & __  “N 

^/p4«>  — <=z,  c<o^s  -4  ^ -y  ^2) 


In  Ref.  1 (pg.192  ) , the  above  egression  has  been  evaluated  for  the  case 
of  rigid  flapping  blades  and  yields, 

r*0-0-* 


-<Z, 


gt-/2 , 
-2 


However,  when  considering  the  flexibility  of  the  blades,  this  expression 
must  be  modified  to  account  for  blade  twist.  Specifically,  this  necessitates 
adding  to  Eq.  (42)  additional  components  of  the  lift,  acting  on  the  blades 
along  the  xQ  axis,  due  to  blade  twist.  These  components  may  be  evaluated 


from  the  expression, 


fa^c*)**/*;  **  fc&J  Or-Yf^  - a,  cost  ^ 

o 

“4  CTOJ-  r 3-/S7  fi/Jc/O 


where 


A&  — — <C  ^ SjO 
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Evaluating  the  various  terms  of  Eq.  (44)  there  results, 
~ Zzr- 


- ~±^  JX*CL,  <& 


£ir  /? 


Ur^r/s? ^b'S'  - -^r  ^ -O^cz,  to  -+-^r^S14;  & 

o o' 

+-±s«  a*  -tl/fA  & -J-s,  a*  (47) 


All  other  terms  equal  zero.  Consequently, 


UHxq) 


= # 

blade 

flexibility 


(48) 


Since, 


frigid  * <aBto) 


blade 

flexibility 
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Adding  Eqs.  (43)  and  (48)  yields, 


+(&,  ■*—<,, )fa. 

-**(*&■*)-  ^ ( 

+ ^ *7-  V^w/ 


where  a = lfJ'bcjQ.R 
2 ^ 

The  equation  for  horizontal  translation  becomes, 

-*-X6r  Mu#  -/■  ey(/  /&%,  V-  o<,  /yL</  -/-<ZZ,  /Yet., 

where, 

Hjux*=  _n .s?7rj 

*\ 

7u#  » y=r  C&J 

Hx  - sP?/?  " 


^ - 


£2  g2>  -/2-  /P^« 
-91- 


/44*  — — 


= 


<X-£Zj*-f2.  Yr^fs.  y.-2 


v 


u.  S^-p-Yi 

St** 


— <^<=e  = 


V 7~ 


s=t  — <ZLdZ£ 
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^ -. 


-S=3g  x,  ^ 

a= »/*P«2 

clo-^  <^\  a.  ./2  S]f 


//'g,  — _ _ Aa, 


Equation  of  fuselage  Pitch 


The  motion  of  the  helicopter  in  pitch  about  the  yQ  axis  may  be  described  as, 


«=*,  ^ 7~/?o</  = <0  (51) 

By  use  of  Eq.  (49);  the  above  equation  may  be  written  in  the  form, 

where 

/? //c*x 

* -Z>« 

+ a AO, 

My**,  -*  ^ ^-/  (52a) 

4*  = ^ ^ 

/^4  - ^4 

^ ^ 
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The  motion  Of  the  disturbed  helicopter,  under  the  assumptions  made,  have 
now  been  fully  specified  by  eight  equations,  two  steady-state  and  six  disturbed 
equations  of  motion.  The  latter  consist  of  two  blade  flapping  equations,  two 
blade,  torsion  equations,  the  equation  of  horizontal  translation  and  the 
equation  of  fuselage  pitch.  However  it  is  immediately  obvious  that  the 
"solution  of  the  six  equations  would  be  extremely  tedious  and  time  consuming. 

It  would  therefore  be  desirable  to  simplify  these  equations,  if  possible, 
without  any  great  loss  in  accuracy.  Such  a simplification  can  be  accomplished 
following  the  procedure  presented  in  Ref.  2. 

Simplified  Equations  of  Motion 

The  quantitative  determination  of  the  response  of  a helicopter  may  be 
greatly  facilitated  by  simplifying  the  blade  equations. 

In  Ref.  2,  it  has  been  shown  that  the  flapping  blade  motions  consist  of  . 

high-frequency  heavily  damped  oscillations.  The  order  of  magnitude  of  these 

oscillations  was  such  as  to  suggest  that  their  influence  on  the  motion  of  the 

helicopter  would  be  small,  i.e.,  the  blade  response  following  a disturbance 

is  rapid  (within  a couple  of  rotor  revolutions)  compared  to  the  long-period 

fuselage  oscillation  which  is  of  principal  concern.  In  fact,  the  disturbance 

hO 

of  the  flapping  motion  is  reduced  to  a tenth  of  its  initial  value  in  ; 

£WT 

revolutions  of  the  blade.  Based  on  this,  little  error  appears  to  be  introduced 
by  neglecting  all  acceleration  terms  proportional  to  a1  and  b . It  is  also 
realized  that  is  approximately  equal  to  , since  approximately 
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equals  o(f  , and  both  quantities  oscillate  at  the  same  frequency.  However, 
since  small  differences  in  and  o<(  introduce  major  contributions  to  the 
response  of  a helicopter  it  cannot  be  assumed  that  a^  equals  o^/  . The  co- 
efficients of  b^  appear  small  and  since  b^  itself  would  be  small,  its 
contribution  may  be  neglected. 

From  the  characteristics  of  typical  rotor  blades  it  may  be  established 
that,  I5«  I1  ; Ig<<  I1  j Ig<<  I1  } I3<<  Iy<<  Ig  » ancl 

I8<<  I2  * 

S’ 

Accepting  the  above  simplifications  appears  not  to  affect  the  long-period 
oscillations  to  any  marked  degree . However,  under  these  simplifications 
the  equations  of  motion  of  the  blades  about  the  flapping  hinge  will  permit 
direct  solution  for  a.^  and  b^  in  terms  of  the  other  degrees  of  freedom. 

By  extension  of  the  above  arguments  to  the  case  of  the  blade  torsion 

99  • •»*  • 

equations,  the  contribution  of  the  dynamic  terms  K jK  3 ^'2.  and 
may  be  recognized  as  small.  Thus,  the  simplified  blade  flapping  and 
torsion  equations  are  as  follows. 

Simplified  Blade  Flapping  Equations 
constant  term  (steady-state), 

x)  (53) 

cosine  term, 

4 - -sU* /%ur  % sty*;  (54) 
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where 


x rxi 


— — -j?~  X° 


- _Q_ 

^yr,=-yZ 


(54a) 


sine  term. 


&/  — *=’*✓  —«*<,  /W/d,  ~^~  ~~ 


where 


M'yux  =■ 

- -& 
*y*  - - 


{£*>  * 


Simplified  Blade  Torsion  Equations 
constant  term  (steady- state), 

y _ 

60  ST* 

cosine  term, 


A~ 


& 


where. 


K- 

72,= 

7a  ■ 


v6^ 


^2^  *4 


£ 

o 

Ca 


/?* 
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(55) 


(55a) 


(56) 


(57) 


(57a) 
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sine  term 


ZdSr. 


7~^ 


7^ 


/* 

A/ 


where 


7~~sC*<--X <77^  T*4  Ao J 

7*'  - ^ 

/cFV  J3" 

7r^y  =-s ?-2^  — /Q_ 

72*  - ^ 


(58a) 


^'=  -^-  -%3^ 


/r3', 

7-4--<ZA. 

It  is  now  possible  to  solve  for  the  quantities  a^,  b^,  ^ and  & in 
terms  of  the. two  degrees  of  freedom,  fuselage  translation,  and  pitch. 

This  results  in  the  following  values, 

£?/  - ***</  - ^ Z7- 

^ •»■«•  -/**.  * /$  — A/ A ^59) 

^ ® /^r  ^ <*£T” 

/a?  — — : A^-xArf •-+  es></  A** 

where,  letting  ^ ..7j£/-.  /ty ^ 

Mi  - ^ ^V/4 
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Substituting  Eqs.  (59)  into  Eqs.  (50)  and  (52),  the  simplified  equations 
of  fuselage  translation  and  pitch  are  thereby  obtained. 

Simplified  Equation  of  Horizontal,  Translation 

Ar  -*l<=A'/  //*(,  V-es^  jg?? ///■'  (60) 

where 

//tlx  ~~  /4  — //<%  £ 

m 

&»,  - ^ 

_ (60a) 

— H£.,-AC,A-+A'+-7+//^  /=> 

M*,  m Mx,  - AC,  =-  7~ 

= /SB,  + AC,  - 7- 


Simplified  Equation  of  Fuselage  Pitch 


% 


(61) 
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where 

M*>mx  ” -Sty,#  £- 

W'S,  “ Af^/ 

- /?  /K*, 

- /> 


(6la) 


Response  of  a Helicopter  to  a Control  Step  Input 

The  response  of  a helicopter  to  a control  input  may  now  be  easily 
determined  by  use  of  Eqs.  (60)  and  (6l)  and  Operational  Calculus . The  general 
procedure  for  the  determination  of  the  fuselage  pitch  response  is  briefly 
discussed  below.  For  the  response  in  translation  (Mr  ),  and  a detailed 
description  of  the  use  of  operational  calculus  for  stability  analyses,  use  of 
Ref.  1 is  suggested. 

Taking  the  Laplace  transform  of  the  above  two  equations  there  is 
obtained, 

7 * - /?/3  ^ 

faj  A&  «*/  ~ /?2-1 
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where, 

*“  rfu*  4 y-  Abx 

^ 4^,  > 

- /42J" 

^■*3  - >4 


The  transfer  functions  is  given  by, 


S*j6\i  » 


A//  Z^/J  j 

■4g/  /4g;ri 


in  which 


where 


4 - 


4,  - 

47- 

4 * 

K - 


A/  Ai 

4?  4 ■/•  V A/  /A  y A0 

'^<r  - /A 

'V**-  /%/  ^ ^5,  /^L,- 

- .Vw,  Mpeji,  w 7“ ^fo^;v 

- ^7 
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(63) 


(64) 


(64a) 
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According  to  Routh ' s stability  criterion,  dynamic  stability  only  occurs 

if  the  coefficients  b , b , b , and  b are  positive  and  if, 

2 2 1 o 

b2  bx  > b^  bQ  (64b) 

For  a typical  helicopter,  b^  is  negative,  and  therefore  is  inherently 
dynamically  unstable.  By  inspection  of  the  terms  involved  in  the  above 
equations,  it  can  be  seen  that  the  major  effect  of  blade  flexibility  and 
unbalance  is  shown  in  the  factors  My^  and  , and 

consequently  bg.  The  values  of  b^,  bg  and  bQ  remain  essentially  un- 
affected. For  mass-overbalanced  or  aerodynamically  under-balanced 
blades  and  H-^  are  increased  and  vice  versa.  This  is,  therefore, 

equivalent  to  increasing  or  decreasing  the  "heaviness"  (coriolis  effect) 
of  the  blades,  which  is  associated  with  damping  in  pitch  of  the  rotor. 

Expanding  Eq.  (63)  there  results 


c/y  A 

<£3  + 

where 

A — //y£~X-  — /rfpOyUje  J 

* i 

} 


(65) 

(65a) 
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Setting  A - 0 yields,  for  conventional  type  helicopters,  one  real  root 
and  one  complex  pair,  say 


^ =A 

A?  * 

Thus,  the  time  history  of  the  pitching  motion  is, 


(66) 


where 


/r e ^ y- /& g £ y-/r3  zfr 


At  »• 


(67) 


By  use  of  the  boundary  conditions,  = -v  =»  C>  at  t = 0,  the 

constants  K_  and  K may  be  evaluated  and  thereby  ^ fel 
<=■  i >=*. 


A similar  approach  may  be  used  to  determine 


AAxdZL) 

S’, 


However, 

g *4, (a) 


since  the  characteristics  of  the  motion  can  be  seen  by  use  of  either 
or  there  is  no  need  to  present,  herein,  the  solution 
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6.  Procedure  and  Application  of  Method 

To  demonstrate  the  procedure  and  application  of  the  method  derived  in 
Section  5>  a typical  helicopter,  equipped  with  aerodynamically  underbalanced 
blades,  will  be  analyzed.  The  aircraft  considered  is  the  same  as  that 
analyzed  in  Ref.  1,  and  has  the  following  specifications, 

W « 5000# 

R * 2k  ft. (the  aerodynamic  effective  radius) 
b = 3 

c = 1.5ft.  (the  effective  chord) 

J2.  = 20.3  rad/sec 
I = 5^0  slug  ft2 
!2  = 0 

W = 102.2#  (the  weight  of  .each  blade) 

^ 2 

IyQ  - 5775  slug  ft  (including  blade  mass  at  hub) 
a «*  5.75  per  radian 
= .00238  slugs -ft3 
h » 6.25  ft. 

f = 0.018 

7 = S&lB?  = 12.65 

-2- 1 

^ TjEsr?'  ■ ■°-o495 

a2  = 62-7 
hx  = -0.082  ft 
K = 23OO  ft.#/rad 

ch/fc  = -2.1, 
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The  step-by-step  procedure  is  as  follows: 

a)  Based  on  the  blade's  mass  and  structural  characteristics  determine 

(r) . This  may  be  approximated  by  the  first  torsional  mode  shape  of  the 

blade.  For  a uniform  bladej  <p  (r)  - sin  . 

2R 

b)  Evaluate  the  mode  shape  integrals  (pg.  ia.  ) 

SQ  = 0.637R 
s1  = 0.405R2 
= 0.296R3 

S3  = 0.231R^ 

c)  Determine  the  steady  state  values  of  ©Q>  A , and  jo  from  the 
equations. 


(53) 


(37) 


(56) 


thus, 

/£,=  0.1^5 

Qc  = 0.213 

jC  - -0.06 
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d.)  Evaluate  the  stability  derivatives  Eqs  (50a)  (52a),  (54a),  (55a), 
(57a),  and  (58a). 


FAr 

* 75,600 

\ 

= -16.3 

= 1958 

\ 

* -4245 

* -970 

Hr, 

= 3769 

= -2175 

H* 

= -1761 

Hax 

= -209.1 

* 2175 

Hal 

* 2825 

My 

°/Vr 

= 12,240 

«y06( 

= -26,530 

= 5775 

My°r, 

= 23,559 

= 17,656 

***<>* 

* -11,004 

Hsrcut. 

- -1307 

» -13,594 

*V, 

= 17,656 

Myo^; 

= -101.9 

^ - 

-6350 

T*x 

= 1008 

Tax  = 

143.9 

= -2921 

Tbi  = 

-2921 

T’a^ 

= 0 

TY,  = 

-6761 

T'a 

1 

= -2921 

TVa 

= -6761 

= 2921 
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= -0.00765 

4/ 

= 0.3272 

= -0.1933 

H&, 

= -0.0623 

= 0.0493 

**  ' 

= -0.924 

- -O.924 

-- 

e)  Determine  blade  motion  influence  factors  Eq.  (59)  • 


L 

- 0.6008 

- 

N.  a 0.6008 

D 

= 0.3154 

A = -0.0127 

F 

= -0.1037 

H a -0.1772 

M 

- -0.0128 

J = 0.0493 

P 

« -0.0449 

C a -0.0055 

E a 0.0174 

Evaluate  stability  derivatives 

Eqs  (60a)  and  (6la) 

Hu* 

X 

= 75,567 

“V,  ■ -208-2 

a 2,009 

=!2,554 

II 

1 

VO 

-J 

O 

= 5775 

55, 

« 371.9 

= 2324 

= -5000 

M^O  a 31,250 

■l 

«7. 

a 5000 

g)  Determine  characteristic  equation  coefficients  Eq,.  (64) 
bl  * 436.2  x 

J 6 

1>  ■ 199-5  x 10 

_ 6 
b « -1.041  X 10 

b"  « 62.TT  x 106 
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j 


,/f 

( 

V.. 1 

t 1 
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h)  Determine  roots  of  characteristic  equation,  A = 0. 

\x=  O.I36  + 0.422  i 
/j*  -0.730 

i)  Evaluate  K, , K , and  K , Eq.  (67) 

A . 2 3 

\ » 5.833 

k2  = -5.833 

k3  = 11.969 


j)  Thus,  the  response  is. 


c*,  (t)  . -0.730t  „ - 0.l36t  _ 

■ B 5*833e  -5*833e  cos  0.422t 


Bn 


0.l36t 

+ 11.969e  sin  0.422t. 


A graph  of  this  response  is  shown  in  Fig.  5. 
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7.  Effect  of  Blade  Flexibility  on  Automatic  Control  System  Requirements 

In  general,  for  the  rapid  subsidence  of  a disturbance  of  a helicopter, 
periodic  control  displacements  (which  tilt  the  rotor  tip-path  plane 
relative  to  the  rotor  shaft)  in  phase  with  fuselage  attitude  («*, ) and 
pitching  velocity  (damping)  are  required  (Ref.  6).  The  first  corresponds 
to  a kind  of  static  stability  and  decreases  the  period  of  the  motion. 

The  latter  increases  the  period,  and  when  applied  alone  tends  to  prevent 
a further  increase  of  the  disturbance.  Moreover,  the  latter  has  a 
remarkable  stabilizing  effect  when  combined  with  the  first.  By  proper 
combination  of  the  two  any  helicopter  can  be  stabilized. 

Since  it  has  been  shown  that  blade  flexibility  and  unbalance 
noticeably  affects  the  damping  in  pitch  of  the  helicopter,  it  is 
reasonable  to  expect  that  these  factors  would  also  affect  the  automatic 
control  system  requirements  for  stability.  To  demonstrate  this  fact, 
consider  a helicopter  equipped  with  an  attitude  gyro  system  so  as  to 
yield  control  displacements  in  phase  with  the  fuselage  attitude  ), 

For  this  case,  Eqs.  (60)  and  (6l)  are  expanded  to, 

/fax  +*Cx  /fax  / )-<,  = 72^  (68) 

^ ■*/  (69) 

^ ^ ~ 

° RESTRICTED 

46. 


RESTRICTED 


The  coefficients  of  the  characteristic  equation  become, 


(70) 


** 

6,  *■  7~?9  Ate*  y-  7~  {/-£#) 


7~  Sf^c^r 


It  is  noticed  that  while  pitch  damping  (MyQ^,  ) affects  the  coefficient  bg, 
automatic  attitude  control  (B*j  affects  the  coefficient  b^.  The  stability 
condition  Eq.  (64b)  may  thereby  be  satisfied.  It  is  interesting  to  notice 
that  without  attitude  control  0*4  = 0),  stability  cannot  be  obtained.  In  other 
words,  increasing  the  pitch  damping  of  a helicopter  can.  only  lead  to  a 
neutrally  stable  craft. 


To  quantitatively  show  how  flexibility  and  unbalance  affect  the 


automatic  control  system  requirements,  assume  a system  in  which  Ba  = O.09. 


The  coefficients  (Eq.  70)  for  various  values  of  mass  unbalanced  blades 


become. 


0 

1.4 

2.2 


436  x ICC 
436  x 105 
436  x 1(> 


b 

■ 2 


106  x 10S? 
393  x Mg 
558  x 10° 


218  x ioi 
218  x 10° 
218  x 106 


b 

o 


62.3  X ic£ 

68.6  x 10° 

71.7  x 10° 
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The  frequencies  and  damping  factors  of  the  fuselage  oscillations  are  therefore, 


i.2  nr 

Damp,  factor 
..  <"i) 

Frequency 
. .(ni>. 

Time  to  half 
amplitude  (sec) 

Period 

(sec) 

K 

0 

0 

0.705 

OO 

8.9 

1.4 

-.175 

0.500 

4 

12.5 

2.2 

-.167 

0.420 

4.2 

15.0 

Thus,  a neutrally  stable  machine  of  short  period  results  when  balanced  or 
rigid  blades  are  assumed.  The  machine  becomes  more  stable  as  the  degree  of 
overbalance  or  flexibility  increases . Finally  a very  satisfactory  stable 
oscillation  of  long  period  is  obtained.  In  the  case  of  underbalanced  blades, 
unstable  oscillations  would  result.  Thus,  consideration  of  blade  flexibility 
and  imbalance  can  often  mean  the  difference  between  a satisfactory  or 
unsatisfactory  automatic  control  system  design. 
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FIG.  7 Coordinate  System 
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